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Abstract. In this work we study the solutions to some fractional higher-order 
equations. Special cases in which time-fractional derivatives take integer values 
are also examined and the explicit solutions are presented. Such solutions can 
be expressed by means of the transition laws of stable subordinators and their 
inverse processes. In particular we establish connections between fractional 
and higher-order equations. 



1. Introduction 

The aim of this work is to investigate the solutions to some boundary value 
problems involving the equation 

(1.1) (D?+D£)u = o, xen ,t>o, Vj e (o,i)un, j = 1,2 

where fta — 17 U {0} and 17 = (0,oo). The symbol D v z stands for the Riemann- 
Liouville fractional operator and will be better characterized further in the text. 
Here the fractional powers uj, j = 1, 2 can take both real and integer values in the 
interval (0, 1) U N. Thus, we study fractional- and higher-order equations. 

Time fractional equations have been investigated by many researchers. The con- 
text in which such equations have been developed is that of anomalous diffusions 



or diffusions on porous media. In the works by Nigmatullin 22] ; Wyss |27| ; Schnei- 



der and Wyss 26]; Kochubei [14J the authors dealt with second-order operators in 



space and gave the solutions in terms of Wright or Fox functions. More recently, 



in the paper by Orsingher and Beghin [24| some interesting results on the explicit 



form of the solutions have been presented. Moreover, in the paper by Meerschaert 



et al. 2.1] and the references therein the reader can find some interesting results on 
abstract Cauchy problems and fractional diffusions on bounded domain. 

In the present paper we deal with higher-order heat type equations in which the 
derivative with respect to t is replaced by the fractional derivative and therefore 
we obtain higher-order fractional equations. Such equations have been studied by 
Beghin Q who has presented the solutions in terms of inverse Fourier transforms. 
In that paper, the author pointed out that such solutions can be viewed as the 
transition laws of compositions involving pseudo-processes and randomly varying 
times. Furthermore, we establish some connection between higher-order equations 
which lead to pseudo-processes (or higher-order diffusions) and, the time-fractional 
counterpart to those equations which lead to either stable or their inverse processes. 
In the matter of pseudo-processes (see Section |4|) we refer to the paper by Krylov 



Date: June 13, 2011. 

2000 Mathematics Subject Classification. Primary . 

Key words and phrases. Fractional equations, Higher-order equations, Subordinators, Inverse 
processes, Pseudo-processes. 

1 



2 



MIRKO D'OVIDIO 



15 1; Hochberg [l3|; Funaki Q; Orsingher (23|; Hochberg and Orsingher 11 1; Lachal 



In sections I2|3I and IH some preliminary results are presented whereas, the 
main results of this work are collected in section [SJ In particular, we obtained the 
solutions to the equation (II . If) in the following cases: 

a) V\ 6 N, 1/2 6(0,1], 

b) v x 6 (0,1], ^eN, 

c) 14 6(0,1], ^2 6(0,1]. 

The special cases 

al) v\ = 1, v 2 = 1/n, neN. 

51) z^i = 1/n, n 6 N, i*j = 1, 

represent the fractional counterpart of the higher-order equations of order n. Sec- 
tion |5] is devoted to the explicit representation of solutions in some particular cases 
involving fundamental equations of mathematical physics. 

2. Inverse processes 

Let ip = f(x,t) denote the distribution of a Levy process Xt, t > on K™ for 
which 

(2.1) Ecxp-^X t =exp 

The Levy process Xt, t > represents the stochastic solution to the equation 

(D 0+it -A)<p = Q, D(A) 

where 

£>(.A) = (/ 6 Lj oc (R n ) : / |/(£)| 2 (1 + **(£)) ^ < oo 

We have used the familiar notation in which /(£) = .F[/ (•)](£) stands for the Fourier 
transform of /. 

A stable subordinator, say S) u (t), t > 0, is a Levy process with non-negative, 
independent and homogeneous increments, see Bertoin |3j and, the cc-Lapace trans- 
form reads as 

(2.2) Ecxp-A-Mt) = C[h y {-,t)](\) =exp-t\ v 

where h v is the density law of S^ v . Straightforward calculations lead to 

. 7W £ 



(2.3) = ier ( -* 2 1; 

According to the literature, we define the process £ u {t), t > as the inverse to the 
stable subordinator Sj v and for which Pr{Sj v (x) > t} = Pr{Z u {t) < x}. Such an 
inverse process has non-negative, non-stationary and non-independent increments 
(see Meerschaert and Scheffler |20|). The law of £ u can be written in terms of the 
Wright function 



z" 



, ^ fc! r ( ak 



as follows 

(2.4) l v (x, t) = i-"^,!-,, (-art - ") , x 6 fi 0> t > 0, v 6 (0, 1). 
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From (|2.4j) we immediately get the Laplace transform 
(2.5) Eexp-A£"(t) = £[/„(•, t)](£) = 

in terms of the Mittag-LefHer function E a (z) — E a ^i(z) where the entire function 



z k 



k>0 

is the generalized Mittag-LefHer for which 

poo \ a ~P 

(2.6) J e - Xz zP- 1 E a .p{-iz a )dz = ^ r ^. 

The function l v is the density law of the inverse process 2, v . The governing equa- 
tions of both processes introduced so far can be written by means of the fractional 
operators (|A.lj) and (|A.2|) as we will show in the next section. In particular, from 
the symbol (|2.3p . we recognize that 9) v is a stable process (positively) totally skewed 
(see e.g. Zolotarev [28]). Thus, as for stable processes, we expect a fractional op- 
erator in space as well. 

We state the following relevant fact. 

Lemma 2.1. The following holds 

(2.7) -Ut,x) = hJx,t), xen,t>0. 
x 

Proof. The (x, t) — Laplace transforms of h v writes 

/•oo y* oo />oo i 

(2.8) / er xt / e-Z x h u {x,t)dxdt= / e^e"^' 'dt 



This is because of the fact that Eexp— £$j u (t) = cxp—t^. From (|2.6j) we obtain 
the Fourier-Mellin integral 

(2.9) C[K(x, -)](A) = ^ J |^ = x v ~ 1 E v , u {—\x v ). 



Let us assume that formula (|2.7[) holds true. Thus, the t-Laplace transform of h v 
can be evaluated as follows 

I" 00 t 
C[h v (x,-)](X) = e- xt -l v {t,x)dt 



x 

1 d r °° 



= - -— / e- xt l v {t,x)dt 
x d\ J Q 

= (by & = -\^E u {-\x») 
=x v ~ l E v ^{—\x v ) 

which coincides with (|2.9[) . We can obtain the same result by considering (|2.2[) and 
the fact that 



roo i-oo \u— 1 

(2.10) / e- £ - x C[l v {x 1 -)]{\)dx = \ v - 1 e-t x exp-x\ v dx = — 

Jo Jo ^ 



The Laplace transform C[l v (x, )\{X) can be easily carried out from (|2.4[) . □ 
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3. Fractional equations on fi 

It is well-known that the law h u of the stable subordinator f) v satisfies the 
fractional equation 

(3.1) (A)+,t + D V Q+X ) h v =0, xefl,t>0 

with initial condition h u (x, 0) = 8{x) whereas, for the inverse process £ v , we have 
that 

(3.2) (D» +t + D 0+tX ) I, =5(x) v ^_ v y xefl o ,t>0 

with initial condition l v {x,Q) = S(x). The fractional operators appearing in (|3.1|) 
and (|3.2[) are the Riemann-Liouville fractional derivatives defined in (|A.1[) and 
(|A.2|) . The initial and boundary conditions for the equation p.lj) can be written as 

h u (x, 0) — S(x), x £ fl, 

K(o,t) = o, t>o, 

whereas, for the equation (I3.2[) . we get 

(3.4) (D» +A + D +,x) l v =0, xeCL,t>0 

subject to the initial and boundary conditions 

l v (x,Q) = S(x), xefl, 
l u (Q,t) = $ 1/ (t), t>0. 

In (|3.5j) we considered the function 

1 



(3-3) 



(3.5) 



r(i-o) + 

where z+ a = H(z) z~ a with a ^ 1, 2, . . . and, H(z) is the Heaviside function. Since 



*a(«) 

® a {z) € we have that 

(3.6) £[$ Q (.)](C)=C Q - 1 . 

From the discussion made so far, the problem of finding solutions for (|3.4I) is to 
trace back through the study of boundary values for l v . This approach will turn out 
to be useful when we study the solutions to higher-order equations which represent 
the higher-order counterparts of both (|3.4[) and (|3.ip . 

The problem of solving (jl.ll) with fractional powers Vj £ (0, 1), j — 1, 2 can be 
approached by introducing the process Srf 1 (Z Vl (i)), t > driven by the law 

(3.7) <)), ieO , t>o, ^e(o,i), i = i,2 

(see e.g. [I; [H; [l9| and the references therein). For v\ = v 2 = v the law (|3.7[) takes 
the form f Vj „(x,t) = t -1 /„(t _1 a;) where 

„ , . 1 x v ~ 1 sin7ri/ „ 

U{x) = - — xen ,t>o 1 1/6(0,1) 

7T 1 + 2x I/ COS 7TI/ + £ 

and (£"(*)) '= t x iSj v (t)/uSj v (t), t > where^C*), j = 1,2 are indepen- 
dent stable subordinators (see for example [3; [H; E3)- We notice that the ratio 
i$) u {t)/ 2^)' / (t) is independent of t. The governing equation of the density (|3.7|) (see 
e.g. [6]) is written as 

(3.8) (£>ft. (t +D^ +tX )j Vl ^{x,t) = S(x) *^ , xeCl Q ,t>0 
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with f VllV3 (dSlo,t) = and f Vl<V2 (x,0) — S(x) or, by considering (|A.4|) . as 
(3.9) (^ + £o+,*W" 2 (M) = 0, xen ,t>0. 

In the next sections we will study the remaining cases in which the powers Uj, 
j = 1,2 can also take integer values. For this reason we will give a short introduction 
on pseudo-processes. 

4. PSEUDO PROCESSES 

According to the literature we define the pseudo-process -Xj , t > 0, n > 2, 
which is a Markov pseudo-process with law satisfying the higher-order heat equation 

(4.1) (J^ - v n = 0, i£l,f>0,«>2 

where n n = (— l) p+1 for n = 2p or n n — ±1 for n = 2p + 1. Such a process is 
termed " pseudo- process" because of the fact that the driving measure is a signed 
measure. We refer to the interesting work by Lachal [l6| and the references therein 
for an exhaustive discussion on this topic. Here, we only recall that, for a given 
order n > 2, the solution to (|4.ip can be expressed in terms of its inverse Fourier 
transform 

1 r + °° 

(4.2) v n (x,t) = — e -iCx+ Kn {-iQ t dC 

when it exists. For n = 3 we obtain a solution in a closed form. Indeed, the 
solution to the third-order heat equation with K3 = ±1 can be written in terms 
of the well-known Airy function (see, for example, Lebedev [l_8] for information on 
this function) and we obtain that 

(4.3) v 3 (x,t) = -^=Ai , xeR,t>0. 

The reader can consult the paper by Orsingher [23| for further details on this case. 
The following result can be viewed as the higher-order counterpart of Lemma 

i2~n 

Lemma 4.1. Let Wj, j = 1,2 be two solutions of (|4.1I) . We have that 

(4.4) wi(x,t) = ~w 2 (x,t), x G M, t > 0. 



Proof. First, for the solution (|4.2j) . we prove that 
(4.5) (2^ + ^1)^=0. 

Let us write 

w„(C,t) = ^[«n (•,«)] (0 = J^ x v n {x,t)dx. 

We have that 

r T 'o+,x 



^Sr>„(-,t)](C) = HC) n - 1 «n(C > *) 



and 



J-[a;«„(-,t)](C) = ~«n(C,*) = -nK^C-iCr- 1 ^^,*). 



6 



MIRKO D'OVIDIO 



By collecting all pieces together we prove the identity (|4.5|) . The next step is to 
evaluate the following derivatives: 

d n w\ n d n ~ 1 W2 x d n W2 
dx n ~ 1 dx n - 1 + 1 dx n 

and 

dw\ x x d n W2 

where we used the fact that W2 solves (|4.1j) . By summing up such derivatives we 
get that 

dw\ d n wi x nd n ~ 1 W2 



dt Kn dx" ~ t* W2 Kn t'dx n 



-1 



d n 1 w 2 x 
t \ nKn ^x^ + l W2 

where W2 solves (|4.ip and thus, the identity (|4.5p is in order. We obtain that 

dw\ d n w\ 

and therefore w± solves the higher-order equation (|4. 1 [) . □ 



Remark 4.2. From Lemma 14. 11 by applying m-times the identity (|4.4j) we imme- 
diately obtain that 

x m 

w(x,t) = —v n (x,t), i£l,f>0,meN 
solves the equation (14. ip . 

5. Fractional higher-order equations on O 

The solutions to the fractional problem involving higher-order operators in space 
of a general order n £ N on the whole real line have been investigated by Beghin 
. In that paper the solutions are presented as the transition functions of pseudo- 
processes with randomly varying time T a - In particular, see [2|, Theorem 2.3], the 
solutions to 

jgzu(x,t) = K n -§-^u(x,t), x £ R, t > 
u(x, 0) = S(x) 



U a (x,t) = / p n (x,s)v 2a (u,t)du 







(5.1) 
coincide with 



where V2 a is the Wright function (|2.4[) . In such problems only the initial conditions 
are required and the solutions are presented in terms of the inverse Fourier transform 
(|4.2[) . Indeed, the functions p n are the solutions to (|4.1j) . 

Our aim is to investigate fractional higher-order equations on the semi-inhnite 
interval Hq with suitable boundary conditions. Due to the fact that /„(0 + ,i) = 
< oo for t > 0, if we are looking for solutions which are symmetric, then 
we can extend these results to the whole real line without effort. It is enough to 
consider the function l„(\x\,t). 

We state the following result. 
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Theorem 5.1. For v € (0, 1], n € N, we have that 

(5.2) (D%_ <x - D» 0+tt )U = 0, xen,t>0 

with conditions (I3.5[) and 

(5-3) D%_J$(x,t) 

Proof. We write the Laplace transforms 



(t), < k < n. 



e~ M U(x,t)dtdx = 



Keeping in mind formulae (|A.3[) and (IA.7I) we evaluate the Laplace transform 



i-iy 



it) 



k=0 



where, once again from (|A.3|) . formula (15.31) has been rewritten as 

D k Q+ Ju{xM_ n+ = {-l) k <b„_ ih±11 {t). 



From p.6p and the linearity of the Laplace transform we get that 



(5.4) 



p oo n — 1 

/ e- At V(-l/£) fc $ Kfc+1) (t)dt 
J ° fc=0 



n-l 



n-1 



fe=0 fe=0 



2- (r - (-i)"A ,y ) 



Thus, from (|A.6j) and ()5.4|) . the equation (|5.2p takes the form 



o = -A^(e,A) + (-iy 



r%(e,A)--^-^(c-(-i) n A*) 



((-in«-A)%(e,A) 



At 



£ + A~ 

-a-irc-A) 



and immediately we get that 



%(£,A) 



At- 



e + A. 



which is in accord with (j2.10j) . This concludes the proof. 
Corollary 5.2. For v — l/n, n E N we have that 
(5.5) (D Q _, t +D%_ jX )l v =Q, xen,t>0 

with conditions (|3.5[) and 



(5.6) 



Dq- x lv(x, t) 



x=0+ 



is the higher- order counterpart of the equation (|3.4 
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Remark 5.3. We notice that $i(t) = and thus, for v = l/n we have that 

D^J u (x,t)\ x=Q+ =0. 

Furthermore, 

D%_ !X l v (x,t)\ x=0+ = B& tt Mt) = D^M0,t). 
We pass to the study of the equations involving the laws of stable subordinators. 
Theorem 5.4. For v € (0, 1], n € N, we have that 

(5.7) (D$_ ti - D% +iX )h*. = 0, xe!l,i>0 
with initial conditions (|3.3[) and 

(5.8) D*h±(x,t) 



t=o+ 



= ^K* + i(a:), < k < n. 



Proof. First, we write 



e- ix h»{x,i) dxdt = h»(£,X). 



By evaluating the x-Laplace transform of the equation (15.71) with the boundary 
condition h v (0,t) — we get that 

(i??_, t -r)M6t) = o. 

By applying standard Laplace technique, in particular from (| A.3|) and (|A.7[) . we 
obtain that 



/ e- xt D% + t h»(x,t)dt = \ n h»(x,\) - \ n ~ 1 5(x) - A™ -1 ^(-1/A) fc $ w= +1 (x). 
J ° k=i 
From (I3.6[) we get 

poo n—l n—1 

/ e^5>l/A) fc ^ +1 (x)dx=£(-^/A) fc 
- 70 fc=i " k=i 

and thus, the (x, i)-Laplace transforms of the equation (|5 . T[) lead to 



0=(-l) n 
-(-1)" 



ra-l 



s (e,A)-A»- 1 -A»- 1 X;(-^/A)* 



fc=i 



n-l 



%(C,A)-A"- 1 ^(-^/A) A 



fc=0 



f%(C,A) 



K-Ar^^A) -r^^A) - (-i)^"- 1 j2(-^/ x ) k 

k=0 

,n-ll-(-^A) n 



=((-A) n -n^(e,A)-(-i) w A f 



l + ^/A 



:((-Ar-nHK,A)- 



A + £i 



Finally, we obtain that 
(5.9) 



%(C,A) = 



A + ^ 



which is in accord with (|2.8[) and this concludes the proof. 



□ 
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From (|5.7|) . for n = 1, we immediately reobtain the equation (|3.1|l . Furthermore, 
a direct consequence of the previous Theorem is the following 

Corollary 5.5. For v = 1/n, n E N, we have that 

(5.10) (Dg_ jt + D _ jSr ) h v = 0, x e ft, t > 
wzi/i conditions (|3.3p and 

(5.11) Dg_ 4 /i„(x,i)| t=0+ = $^+i(x), n>fc>0. 



is ifte higher-order counterpart of the equation (|3.1|) . 

Remark 5.6. The functions Z w and ft.!,, for v = 1/n, n e N can be written as Mellin 
convolution of generalized gamma functions or, for n € 2N + 1, in terms of Mellin 
convolution of the Modified Bessel function K a (see 0;@). 



6. Higher-order equations: further directions 
Let v(-,t) be in C(f2 ) as a function of t. We have that 

(v(x,-), D 0+ ,./i,i (■,*)) =Nn(x,t)-{h±(;t), Do+,.v(x,-)) 

where 

N n {x, t) — v(x, s) hi (s, t) 
For n > 1, the conditions (|3.3p hold true and then we get that Af n (x, t) = and 

(«(«,•), D 0+ ,.hi_{;t)) = {hx(;t), D -,.v(x,.)). 
Thus, the stochastic solution to the equation 

(6.1) (D%_ t + D m Ju m ,„ = 0, ie!l,(>0, rn, n > 1 

is given by the process £™(f}i(i)), t > with law 

U„ ltn (x,t) = (lx(x,-), h±(-,t)). 

From the fact that l v (x, s) ^> 6(x — s) and h u (s,t) '—^ 6(s — t) we obtain that 
u m .i = l\/ m and Ui t n — h%/ n . As a direct consequence we also obtain that 

lim Af n (x, t) = tp(x,t), xGflo,t>0 

n— >1 

which must be understood in the sense of distribution. 

6.1. m = 2. The case m = 2 leads to the process \B(fji t > where \B\ is the 
reflecting Brownian motion. From the subordination principle we have that the law 
of \B(fy~)\ coincides with the folded law of a 2/n-stable process with characteristic 
function flU) and #(£) = \£\ 2/n . Thus, for n = 2 (and rn = 2), we obtain the 
folded Cauchy process driven by the Laplace equation. 
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6.2. 771 = 3. For m = 3 we focus on the following few cases: 

(D 0+ , t + D 3 0+ J Ul = 0, 

(Dl+, t - D 3 0+ Ju 2 = 0, 

(D 3 0+>t + D 3 + Ju, = 0. 

The explicit solutions to the equations above are listed below: 

1 fx ( 2 a; 3 / 2 
u 1 (x,t) = - ] l- t K i 

(which is the solution of (|4.1[) on the positive real line) where (see Lebedev 18L 
formula 5.7.2]) 

= vI-M - I v {z) 
2 sin vk 

is the modihed Bessel function of the second kind (or Macdonald's function) and 

{z/2f k + v 



fc=0 



fc!r(fc + ^ + i) 



is the modified Bessel function of the first kind (see [HI, formula 5.7.1]); 



1/3 / x 3 \ , . . f2x 



,3 



U2( "'' ) = 4V^r eXP l^J W --n:!^ 
where (see Lebedev [TU, formula 9.13.16]) 

W a ,f)(z) = Z P +1 / 2 e~^U(l/2 -a + j3,2/3 + l;z) 

is the Whittaker function (U is the confluent hypergeometric function of the second 
kind); 

2t x-t 
U3[X,t) ~ 3 3 / 2 tt x 3 - i 3 ' 
The solution w 2 can be obtained by considering the integral 

, s f°° , s te ~^ , 

U2\x,t) = / ui(x,s) as 
Jo v Airs 3 

and the formula 6.631 of Gradshteyn and Ryzhik [9j. The solution U3 comes out 
from the integral 



f°° 1 < 3 / 2 / 2 t 3 / 2 \ 



ds 



and the fact that 

■ n{yz)-"(y 2v - z 2v ) 
2 smirv (y 2 — z 

(see 0, formula 6.521]). We observe that 

2 t 



3 K v {ys)K v {zs)ds= ^ZJ1.2_ 5T » + > 0, < 1 



u 3 (x,rj) 



3 3 / 2 tt a; 2 + + < 2 
which agrees, in some sense, with ()3.7j) for ^1 =7/2- 
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6.3. m = 4. For n = 1, the solution to (|6.1|) is the law of the folded iterated 
Brownian motion |B(|.B(t)|)|, t > 0. This is because of the fact that Z1/4 = ii/2°?i/2 
is the density law of £V 2 (£V 2 (t)), t > 0. The solution y(x,t) = u 4 ,i(|a;|, t), x G K, 
t > has been thoroughly studied (see e.g. Allouba and Zheng [l[) and solves the 
generalized equation 

where (§", <p) = 0"(O) for <j) G C C °°(R). In this case we have that 

For n = 2 we have that r(x,t) = u^ t 2(\x\,t) solves the equation of vibrations of 
rods on K and the corresponding process is B(\Si/4(t)\) where B is the standard 
Brownian motion and S a is the Levy process with ^siQ = and a £ (0, 2] (that 
is, an a-stable symmetric process). 

Appendix A. Fractional and Higher-order derivatives 

We recall the Riemann-Liouville fractional derivatives 

1 d n f z 

(A.l) Dg + J(z) = T[n _ a) 1 ^ I (* - s) n - a - l f(s) ds, z G R + 
and 

1 d n f°° 

(A.2) ^ m = -—-—J a (s-zr-^f(s) d s, zeR + 

which hold for n — 1 < a < n, n £ N. Formulae (|A.1[) and (IA.2I) are usually 
refered to the right and left Riemann-Liouville fractional derivatives and become 
an ordinary derivative for a = n G N. Indeed, 

(A.3) D% +<z = (-1)"^% = d n /dz n . 

For the sake of simplicity we will write Dq+. z instead of Dq + z . Furthermore, we 
recall that 



(A.4) D^J(z) = ^l( z ) + J2D k 0+ J(z] 

k=0 



z k- 



T(k-a + lY 

z—0+ 



n — 1 < a < n 



(see Gorenflo and Mainardi [8]; Kilbas et al. [Hi; Samko et al. [25j]) where 
(A.5) _/(*) = _ / _«)»-«-! JLi( fl ) da n-Ka<n 

is the Dzerbayshan-Caputo fractional derivatives. The Laplace transform of (|A.4j) 
is written as 

(A.6) C[D- + J(-)](C)=C a C[f(-)}(0- 

For the sake of completeness we also recall the Laplace transform 



(A.7) £[dz +! J(-)](0 = C"£[/OK0 - E C^D^Jiz) 



k=0 



z=0+ 
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